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Abstract 

Naturally occurring ices lie on both interstellar dust grains and on celestial 
objects, such as those in the outer solar system. These ices are continu¬ 
ously subjected to irradiation by ions from the solar wind and/or cosmic 
rays, which modify their surfaces. As a result, new molecular species may 
form which can be sputtered off into space or planetary atmospheres. We 
determined the experimental values of sputtering yields for irradiation of 
oxygen ice at 10 K by singly (He+, C+, N+, 0+ and Ar+) and doubly (C^+, 
and charged ions with 4 keV kinetic energy. In these laboratory 
experiments, oxygen ice was deposited and irradiated by ions in an ultra 
high vacuum chamber at low temperature to simulate the environment of 
space. The number of molecules removed by sputtering was observed by 
measurement of the ice thickness using laser interferometry. Preliminary 
mass spectra were taken of sputtered species and of molecules formed in the 
ice by temperature programmed desorption (TPD). We hnd that the exper¬ 
imental sputtering yields increase approximately linearly with the projectile 
ion mass (or momentum squared) for all ions studied. No difference was 
found between the sputtering yield for singly and doubly charged ions of the 
same atom within the experimental uncertainty, as expected for a process 
dominated by momentum transfer. The experimental sputter yields are in 
good agreement with values calculated using a theoretical model except in 
the case of oxygen ions. Preliminary studies have shown molecular oxygen 
as the dominant species sputtered and TPD measurements indicate ozone 
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formation. 
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1. Introduction 

Studies of ion-induced processing of astrophysical ice analogues in the 
laboratory are relevant to a variety of different environments, such as the icy 
mantles of dust grains in the interstellar medium and protoplanetary discs, 
together with the surfaces of Solar system objects such as comets. Centaurs, 
Kuiper belt objects, and the icy satellites of the outer planets. Sputtering 
by ion impact can give rise to changes in the chemical, physical and optical 
properties of the ice as a result of both elastic and inelastic collisions. In the 
sputtering process, material is ejected from the surface which may assist in 
the development of thin atmospheres such as those found around Pluto and 
the icy moons of Jupiter and Saturn. An extreme example is Europa, where 
O 2 on the trailing hemisphere of Jupiter’s moon Europa is sputtered by low- 
energy ions to become the dominant component of the moon’s atmosphere 
and neutral gas torus [DEIE], together with causing the hemispherical colour 
variations. 

Sputtering also plays a role in the chemical alteration of distant comet 
surfaces, and contributes to the formation of the hrst coma as these objects 
approach the sun. Recently, the signature of this process was recorded by 
the ROSINA mass spectrometers on-board the Rosetta spacecraft at comet 
67P/Churyumov-Gerasimenko. Small amounts of refractory elements such 
as Na and Si are believed to have been generated by Solar wind sputtering 
of dust grains on the surface of the nucleus [1]. 

O 2 ice is not a major component of the observed astrophysical ices in¬ 
ventory but O 2 ice processing by the slow component of the Solar wind and 
cosmic rays plays an important role in the formation of other species such as 
ozone [5], water and various carbon-oxides |6]. Whilst there have been many 
sputtering studies using water ice as a target; (see [7j and references within) 
sputtering of oxygen ices has seen fewer investigations. Given that cosmic 
ray impacts onto molecular oxygen on dust grain surfaces is one possible 
pathway in the formation of water (e.g., m), the associated loss pathway of 
sputtering warrants further investigation. Depending on the energy of the in¬ 
coming projectile ion, experimental sputtering can be divided into two main 
categories: sputtering by low energy ions (a few keV) and sputtering by high 
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energy ions (>30 keV, up to ~3 MeV). At energies of a few keV ions primary 
lose energy in ‘nuclear’ ballistic collisions, where target molecules are excited 
vibrationally, rotationally and translationally. Above 30 keV energy is lost 
primarily in electronic excitation of target molecules by, e.g., ionisation. In 
between these two energy ranges neither of these two energy loss mechanisms 
dominates. 

Sputtering of oxygen ice by low energy (4-10 keV) H+, H^, and Hg and 
high energy ions He^ and (up to 3.5 MeV) has been intensively studied 
by Ellegaard et ah IHli, and references within. They concluded that oxygen 
ice sputters more efficiently than nitrogen ice by a factor of almost 2, and 
that the sputter yields are proportional to the electronic stopping power, i.e. 
the energy lost by the projectile ion due to inelastic collisions with bound 
electrons as it penetrates the ice. The effect of multiply-charged ions on ion- 
induced sputtering of O 2 has been neglected so far in the literature, with the 
only study at low energy for oxygen ices by Gibbs et al. [10] on oxygen ice 
irradiated at MeV energies by He^+ and He’*' ions. They concluded that the 
yields are a factor of 1.2 higher for He^"*" ions. 

The present work is a laboratory study of sputtering of oxygen ice by 
4 keV singly and doubly charged C, N, and O ions and singly charged He 
and Ar. There have been only a few laboratory experiments on sputtering 
by low energy non-noble gas ionic species such as C, N and O, which have 
the potential to be incorporated in new molecules in surfaces and even fewer 
experiments for sputtering by multiply charged ions. 

2. Experimental details 

The experimental arrangement used in the present work is shown schemat¬ 
ically in Figure An Ultra High Vacuum (UHV) chamber was coupled to a 
low energy ion accelerator m equipped with a 10 GHz Electron Gyclotron 
Resonance (EGR) ion source [T2|. Singly and doubly charged ions of ^^G, 
i4n, 16 q singly charged ions of ^He and ^°Ar, at 4 keV, were focused 
and directed into the experimental chamber (which operated at a base pres¬ 
sure ~ 1 X 10“® mbar). Ion beams were collimated by three 4 mm diameter 
apertures separated by 50, 38 and 26 mm respectively from the substrate 
providing an ion beam of 4 mm in diameter with less than 1% increase in 
diameter for the maximum current density used at the ice surface. A potas¬ 
sium bromide (KBr) substrate, 20 mm in diameter and 2 mm thick, was 
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Figure 1: Schematic diagram of the experimental arrangement with 1) ion irradiation part, 
2) KBr substrate, 3) radiation shield, 4) vapour deposition nozzle , 5) and 6) thickness 
monitor comprising (5) a 405 nm diode laser and (6) a photodiode detector and 7) line of 
sight quadrupole mass spectrometer. 


clamped to an earthed, temperature controlled sample holder (A.S. Scien- 
tihc) surrounded by a gold plated copper radiation shield. 

The substrate temperature was controlled to within 0.1 K at temperatures 
from 8 to 300 K using a calibrated silicon diode. A slot in the radiation 
shield 15 mm high and with a total angular width of 120 degrees enables 
ion beam access to the deposited ices; the slot gives access of ±60 degrees 
with respect to the normal of the substrate surface in the radiation shield. 
The radiation shield temperature was measured throughout the experiment 
by a diode, and was found to be constant at a value of 48 K during the 
irradiation and heating procedure. The cold-head was rotatable through 
360 degrees via a differentially pumped rotary stage. The substrate could 
be rotated to any desired position during the experiment. Additionally, as 
shown schematically in Figure a quadrupole mass spectrometer (QMS; 
model HAL-201 from Hiden Analytical) could be used to monitor species 
emanating from the sample either from sputtering or thermal desorption. 

The ice films were deposited at a pressure of 1.6 x 10“^ mbar with the 
thickness measured by laser interferometry using a 405 nm diode laser and a 
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Figure 2: Reflection and refraction model for an oxygen ice film deposition on a KBr 
window. 

photodiode detector. During deposition, the cold-head assembly was rotated 
such that the vapour deposition nozzle was directed normal to the substrate 
surface, at a distance of 25 mm, with the laser and detector at 45 degrees 
with respect to the substrate normal. Oxygen ice dims of up to 258 nm were 
deposited onto the substrate which was held at 10 K. 

Figure [^shows a schematic diagram of the reflection and refraction for the 
oxygen ice him deposited on a KBr window. 6q is 45 degrees and represents 
the angle of incidence of the laser light with respect to the surface normal, 
and 6i is the angle of refraction at the vacuum ice interface calculated from 
Snells Law Q. riQ is the vacuum refractive index and rii the oxygen ice 
refractive index. 


( 1 ) 


riQ sin 6o = ni sin 9i. 


The intensities Ji to J 5 can be calculated from the rehection and refraction 
laws. The model ht equation is expressed as 



where PS is the photodiode signal in volts, Ji is the rehection intensity from 
vacuum to the ice and J 5 the refracted intensity after the light passes through 
the ice and is rehected back. In the model a value 1.285 for the refractive 
index of the ice gives the best ht. C is a constant related to the photodi¬ 
ode efficiency and T is the transmission coefficient of the ice at 405 nm. In 
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equation (|^, m is an integer, where odd values of m correspond to max¬ 
ima (constructive interference) and even values of m correspond to minima 
(destructive interference) at intervals of 86 nm for this wavelength. 


d = cos 9i 


( 3 ) 


The ice density was calculated using the Lorentz-Lorenz relation na 


L X p 


— 1 
n2 + 2 


( 4 ) 


where L is the Lorentz-Lorenz factor, n is the refractive index of the ice and 
p the density of the ice. A value of L = 0.1294 cm^ g“^ has been used in 
the present work [13] and we assume that L is constant in the visible wave 
range. Using equation (|^ an oxygen ice density of p = 1.38 ± 0.20 g cm“^ 
was obtained. 

Fulvio et ah [ 13 ] measured an oxygen ice refractive index of n = 1.322 
by laser interferometry and used this value to obtain an oxygen ice density 
of 1.54 g cm“^ at 16 K, again using a Lorentz-Lorenz value of L = 0.1294 
cm3 g-i et ah na found a refractive index of 1.25 for oxygen ice at 

20 K. Therefore our measurements are in line with previous determinations 
of refractive index and density. 

After a 258 nm him of O 2 ice was deposited, the cold-head was rotated 
to enable irradiation by the ion beam, at 45 degrees with respect to the 
normal to the ice him surface. The ion beam current was measured by a 
Faraday cup inserted immediately in front of the ice him and the current 
was monitored during irradiation of the ice sample by a 90% transmission 
metal mesh. Both the Faraday cup and the mesh were positively biased at 
50 V to suppress secondary electron emission. The ice him thickness was 
measured by the laser interferometer after a specihc ion irradiation dose and 
this process repeated for further ion doses. 

Figure [^shows the interference fringes from the photodiode together with 
a ht from the model as described above. Measurements of the reduction of the 
oxygen ice him thickness by the laser interferometer enabled the sputtering 
yield Ys (molecules / incident ion) to be calculated. The ion beam was 
calculated to be an ellipsoid of area 17.7 x 10“^ cm^ giving an ion beam 
current density of ~ 0.15 pA cm“2 

Figure shows the photodiode detector voltage as a function of ion dose 
for irradiation of a 258 nm oxygen ice him, at 10 K by 4 keV N"*" ions. Also 
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Figure 3: Laser interferometer photodiode voltage (points) for experimental values and 
model fit (solid line) vs. time during the formation of oxygen ice at 10 K from oxygen gas 
at a pressure of 1.6 x 10“^ mbar. 

shown in Fignre is a model fit (solid line) nsing the valnes of refractive 
index and density determined from the deposition prohle shown in Fignre 
which allows the determination of the ice thickness removed (top axis). 

3. Results and discussions 

Figure 1^ shows the number of oxygen ice molecules eroded at 10 K, as a 
function of the dose of and ‘^°Ar+ ions, all at 4 keV 

kinetic energy. Figurej^is the same, but for doubly charged and 

i6o2 + ions. Note that care was taken to correct the ion acceleration voltage 
for the charge state of the ion to maintain a constant 4 keV kinetic energy; 
acceleration of a doubly charged ion through a potential difference of 2 kV 
gives an increase of 4 keV in kinetic energy whereas a potential difference of 
4 kV is required for singly charged ions. It is clear from both Figures and 
l^that the number of molecules removed varies approximately linearly with 
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Figure 4: Laser interferometric data (points) and sputtering model (solid line) for 4 keV 
N+ on a 258.4 nm oxygen ice film deposited and irradiated at 10 K. 

the ion dose in all cases.A linear sputter rate was observed for all studies re¬ 
ported here. The persistence of this linear relationship shows that for all ion 
exposures there remained sufficient O 2 ice to prevent interactions between 
the ion beam and the KBr substrate. Furthermore, the linearity implies that 
there was no charging effect of the target ice by the ion beam, which would 
have reduced the intensity of the impinging ions. The experimental error 
in the number of molecules removed at each data point in Figures and 
is estimated to be 16% from htting the experimental measurements of light 
intensity to the model of the light interference patterns. The gradients of 
the straight line linear regression hts correspond to sputtering yields, Ys; 
the number of molecules ejected per incident projectile ion. Experimental 
sputtering yields resulting from these data are shown in Table 1 with uncer¬ 
tainties reflecting only the random error in determining the gradients of best 
£t lines to the data in Figures and As mentioned above, in addition to 
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Figure 5: Number of O 2 molecules removed by 4 keV singly charged ions on oxygen ice at 
10 K as a function of ion dose, where upward triangles correspond to ^He+, circles corre¬ 
spond to squares correspond to diamonds correspond to and downward 

triangles correspond to "^°Ar+. 

these random errors there is a potential systematic error of 16% associated 
with the determination of the number of molecules removed from the optical 
measurements. These systematic errors are due to uncertainty in the exact 
refractive index of the oxygen ice and are, therefore, the same for all the 
experiments. Thus, experimental data can be reliably compared from one 
ion to another by taking into consideration only the random errors shown in 
Table 1, but the additional systematic error should be included when the ab¬ 
solute values for sputtering yields are used in models. In tests with variable 
ion beam intensity the sputtering yields were found to be independent of the 
ion fluence. 

The experimental sputter yields have been compared with values pre¬ 
dicted by a model developed by Fama et ah [7], which are also shown in 
Table 1. Fama’s model was originally used to predict the sputtering yield. 
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Figure 6: Number of O 2 molecules removed by 4 keV singly charged ions on oxygen ice at 
10 K as a function of ion dose, where circles correspond to squares correspond to 

i4]\^2+ diamonds correspond to ^^0^“*". 

F(£', mi, Zi, 0, T), following ion impact on water ice at temperatures, T, 
from 10 K to 140 K with ion energies, E, up to 100 keV and projectile ions 
with mass mi, atomic number Zi and incident angle 6] their expression is 

r (E, rn„Z„0,T) = ± + r, x (l + F .xp cos-/ 0. 

( 8 ) 

Here, Uq is the sublimation energy per molecule. A value of 0.07 eV/molecule 
was used in the current work, which is close to the value of 0.09 eV/molecule 
used by Ellegaard et al. mu- This value of 0.07 eV/molecule was chosen to 
maximise the agreement between sputter yields predicted by the model and 
the experimental values. It has been reported that the exact value of Uq de¬ 
pends on the deposition conditions and ice thickness [U] , which may account 
for the difference between the present value of Uo and the value of Ellegard. 
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Table 1: Experimental and theoretical values for sputtering yield of oxygen ice at 10 K by 
4 keV singly and doubly charged ions. The experimental errors shown are random error 
in fitting the data shown in Figures and There may be an additional systematic error 
estimated to be ±16% (see text). 


n p (g cm ls(02) Expt. 15 ( 02 ) Theory 


He+ 

1.285 

1.378 

25 ± 1 

25 

C+ 


5? 

89 ±5 

94 

N+ 


5? 

119 ±7 

112 

0+ 



170 ± 10 

124 

Ar+ 



252 ± 15 

243 

C^+ 



101 ±6 


N2+ 



106 ±6 


02+ 



176 ± 11 



Sn is the nuclear stopping cross-section and Se is the electronic stopping 
cross-section; these stopping ranges have been calculated with the Stopping 
Range of Ions in Matter (SRIM[^ Programme. In this work, the energy in¬ 
dependent formulae giving values of p (a function of the atomic number, Zi, 
of the projectile), / (the order of the angular dependence, function of the 
projectile atomic mass, mi) and a (function of the mass ratio of target and 
projectile, m 2 /mi) have been calculated using the empirical formulae given 
by Fama et al. [7], which are not reproduced here. The formula used by Fama 
et al. [7] to calculate the value of rj was developed specihcally for water ice; 
it is based on a double exponential £t to a number of previously measured 
values taken from Sigmund [12]. We experimented with the £t, which can 
be done for any particular mass ratio of target ice molecule to incoming ion, 
and found very small differences (<8%) between using oxygen ice and water 
ice, so we decided to retain the original formalism by Fama et al. [7]. Note 
that we use atomic mass 16 rather than molecular mass 18 in the model as 
we assume that the primary interaction is between the projectile ions and 
the individual oxygen atoms in the ice. The projectile ion energy, 4 keV, is 
nearly three orders of magnitude greater than the bond dissociation energy 
of molecular oxygen, ~5 eV. Thus, although collisions will be between the 
projectile ion and O 2 molecules of mass 32 the bond between the two oxy¬ 
gen atoms is not strong enough to hold the molecule together in collisions 
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Figure 7: Sputtering yield per incoming ion as a function of ion mass. Solid symbols 
correspond to singly ionised species (^He+ = upward triangle, = circles, = 

squares, = diamonds, = downward triangles) and empty symbols correspond 

to doubly ionised species = circles, = squares, ^^0'^+ = diamonds.) The 

dashed line indicates the theoretical yield as a function of projectile mass, and the solid 
line shows a simple linear fit (slope 6.79) to all experimental data except the oxygen ion 
values. 


where one of the oxygen atoms is struck by the projectile. Furthermore, 
the stopping ranges calculated using SRIM also consider atomic rather than 
molecular masses. Cq is a constant, which is calculated from the differen¬ 
tial cross section for elastic scattering in the binary collision approximation 
[7]; here we have determined a value of 1.808 from Sigmund [16]. In the 
present work the temperature, T, was constant at 10 K, the energy, E, was 
4 keV and the angle of incidence, 9, was 45 degrees. The term ^ exp 
is neglected here as our experiment was carried out at constant temperature 
and energy (for details, see [7|.) 

Figure shows our experimental O 2 sputtering yields plotted against in¬ 
coming ion mass. Overplotted are the theoretical prediction using the model 
just described, and a simple linear £t with a slope 6.79 molecules per unit ion 


12 




Figure 8 : Sputtered species due to irradiation of O 2 ice by Ar+ ions. Lines indicate 
pressure in the chamber due to different species labeled according to their mass (mass 16 
= O, mass 32 = O 2 , mass 40 = Ar, mass 48 = O 3 ). The dominant oscillation corresponds 
to the ion beam being turned on and off. 

mass to the experimental data. The sputtering yield is therefore proportional 
to mass of the incoming ion or the square of the momentum. As can be seen 
in Table and Figure the experimental and theoretically predicted sput¬ 
tering yields agree with each other remarkably well, particularly considering 
the model was developed for water ice. The only exception is the case of oxy¬ 
gen ion sputtering, where we measure a ~45% higher sputtering yield than 
predicted by theory. We note that the effect is seen in different experimental 
runs using both singly and doubly charged oxygen ions, which give similar 
sputtering yields within the experimental error. It is unclear why oxygen ion 
sputtering yields should be higher than predicted by the model; the matter 
deserves further investigation. As this effect is observed only in the case 
where the target and projectile species are the same, it would be interesting 
to see whether a similar enhancement is observed in N 2 ice irradiated by N 
ions. 
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Figure 9: Temperature programmed desorption of O 2 ice following irradiation by C“'" ions. 
Lines indicate pressure in the chamber due to different species labeled by their mass (mass 
13 = mass 16 = O, mass 29 = ^^CO, mass 32 = O 2 , mass 45 = ^^C02, mass 48 = 
O3). 


As can be seen in Figure the sputtering yields increase as the ion 
momentum increases for both singly and doubly charged ions, showing that 
momentum transfer alone is the dominant factor for almost all the species 
studied. In fact, with the exception of the oxygen ions discussed above, a 
linear relation (6.79 O 2 molecules removed per incoming ion per unit mass 
of incoming ion; solid line in Figure to our experimental data provides a 
useful approximation to the more detailed model in the range of ion masses 
considered. The dashed line in Figure [^connects sputtering yields calculated 
with the theoretical model for the highest abundance isotopes of atoms with 
atomic numbers from 1 to 18; H to Ar. Electronic stopping dominates in 
the case of He’*' whereas nuclear stopping becomes more important for the 
heavier ions. Our result can be compared with Ellegaard et ah |H| who studied 
sputtering of oxygen ice by H+, and with energies from 1-5 keV. We 
measure a sputtering yield of 25 molecules per incident He"*" ion (mass 4) 
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which is consistent with their yield of ~30 molecules per 3 keV ion (mass 
3) at higher temperatures (15-20 K) and normal incidence. 

During irradiation, we monitored the sputtered material using the QMS. 

For example, Figure shows the sputtered species observed when O 2 ice 
was irradiated by Ar"*" ions. In all cases, we found that O 2 is the dominant 
ejected species, in good agreement with Ellegaard et ah Other species 

monitored were O (mass 16), O3, (mass 48), and the projectile ion mass. The 
measured ratio of O 2 to O is the same during sputtering as when O 2 is fed 
into the chamber with no beam. This indicates that the O signal arises from 
the dissociation of O 2 in the ioniser of the QMS. Also, the O 3 signal with 
no beam is the minimum detectable pressure in the chamber so it represents 
the zero level. 

In addition, we performed thermal programmed desorption (TPD) of the 
O 2 ice after irradiation with ions, which consisted of QMS monitoring 
of the molecular species desorbed from the ice as the temperature increased 
(Figure]^. We selected to track masses that might be expected to form based 
on the projectile ion and target ice. The TPD showed no formation of CO 
or CO2, but we found an O3 peak at temperature ~56 K which is consistent 
with predictions by Ellegaard et ah [9] and the experimental measurements 
of Ennis and Kaiser [5] . These authors conducted a comprehensive study of 
O 3 formation in O 2 ice at 12 K irradiated by 5 keV singly charged ions of C, 

N and O. They observed that the O 3 yield did not depend on the mass of 
the projectile ions and concluded that electronic stopping regimes dominate. 

In a study of O 2 ice irradiated by 100 keV protons, Fama et ah PH find that 
the formation of O 3 quenches the sputtering yield of O 2 . Our experiments 
never enter this regime, as evidenced by the linear relation between sputter 
yield and ion dose (see Figures and |^. 

4. Conclusions 

We have measured the sputtering yield from oxygen ice deposited and 
irradiated at 10 K by 4 keV singly charged He’*', C’*', N"*", O’*" and Ar"*" and 
doubly charged and Our findings are as follows: 

• The sputtering yields increase linearly with incident ion mass or with 
ion momentum squared for all studied species apart from O. The yield 
varies from 25±1 molecules/ion for He"*" (mass 4) to 252 ± 15 molecules/ion 
for Ar"*" (mass 40) at a rate of 6.79 O 2 ice molecules removed per in¬ 
coming ion per unit mass of the incident ion. 
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• When using C, N and O ions, the experimental sputtering yields are 
the same for singly charged and doubly charged ion impact within the 
experimental errors, as these possess effectively the same momentum. 

• The experimental yields for singly charged ions have been compared 
with a theoretical model adapted from Fama et ah and the results 
are in good agreement within our experimental errors except in the 
case of oxygen ions. 

• Preliminary analysis of the ejected species during sputtering with a 
quadrupole mass spectrometer showed O 2 to be the dominant sputtered 
species. Temperature programmed desorption of the oxygen ice after 
irradiation with showed the formation of ozone, which is in good 
agreement with the findings of Ennis and Kaiser [5] . 
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